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EFFECT OF SUPPORT FLEXIBILITY AM) DAMPING ON THE DYNAMIC RESPONSE 
OF A SINGLE MASS FLEXIBLE ROTOR IN ELASTIC BEARINGS 
Ly R. G. KIRK 1 and E. J. GUNTER 2 


I NTRODUCT I ON 

The study of rotor dynamics has in rscent years, become of increasing 
importance in the engineering design of p^wer systems. With the increase 
in performance requirements of high-speed rotating machinery in various 
fields such as gas turbines, process equipment, auxiliary power machinery 
and space applications, the engineer is faced with the problem of designing 
a unit capable of smooth operation under various conditions of speed and 
load. 

In many of these applications the design operating speed is o~ten 
well beyond the rotor first critical speed, and under these circumstances 
the problem of insuring that the turbomachine will perform with a stable 
low-level amplitude of vibration is extremely difficult. 

At the turn of the century H. H. Jeffcott (I) developed the fundamen- 
tals of the dynamic response of the damped single mass unbalanced rotor 
on a massless elastic shaft mounted on rigid bearing supports. The Jeffcott 
analysis of the single mass model showed that operating speeds above the 
first critical speed were possible and that a low level of vibration would 
be attained once the rotor had exceeded the first critical speed. 

As various compressor and turbine manufacturers adapted the flexible 
rotor design concept in which the rotors were designed to operate above the 
first critical speed, various units developed severe operating difficulties 
which could not be explained by the elementary Jeffcott model. 

Under certain conditions of high spesd operation above the first 
critical speed, such influences as internal rotor friction (2) hydro- 
dynamic bearing and seal forces (3) and aerodynamic cross coupling (4) 
can lead to a destructive nonsynchronous precessive whirl motion being 
developed in the rotor system. 
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B. L. Newkirk and Kimball (5), in their early investigations of 
instabilities in compressors due to internal friction, were able to deter- 
mine experimental ly that the introduction of a flexible support system 
could greatly extend the rotor stability threshold speed. D. M. Smith 
(6) in 1933 was the first to verify Newkirk’s findings theoretically by 
expanding the Jeffcott model with internal damping to include a massless 
damped flexible support system. Recent investigators such as Lund (7), 
Tondl (8), Dimentberg (9) and others (10) have shown that flexible damped 
supports may improve the stability characteristics of high speed rotors. 

The present analysis was undertaken to primarily determine the in- 
fluence of flexible supports on the synchronous unbalance response of the 
single mass Jeffcott rotor, and to optimize the support system character- 
istics so as to minimize the rotor amplitude and forces transmitted over 
a given speed range. The problem of bearing forces transmitted has been 
examined by various researchers, (II, 12, 13, 14) they have shown that a 
significant reduction in the forces transmitted can be achieved by the 
proper design of the bearing support system. 

For example, Den Hartog (15) has shown that the tuned vibration ab- 
sorber will greatly reduce the response of the forced vibrations of the 
two-mass system. The following analysis parallels this approach for the 
case of a single mass rotor excited by an unbalance load. The analysis 
presents an analytic study of the tuned damper support system similar to 
that employed by Broch (16) and also presents a generalized study performed 
on the digital computer to obtain optimum support damping to produce the 
best response of the rotor over a wide speed range. It is well known that 
a damper support system can improve the vibration characteristics of a 
rotating shaft and various investigators have considered the problem 
either from the standpoint of a continuous elastic system or as a series 
of lumped masses (17-23). 

Although the results presented in this paper apply specifically to 
the single mass Jeffcott model, the optimization procedure may be readily 
extended to more complex multi-mass rotor bearing systems by employing a 
finite element rotor digital computer program similar to the procedure 
presented by Lund in Ref. 24 or by using the procedure as outlined In the 
paper presented by Crook and Grantham (25) on the vibration analysis of 
turbine generators on damped flexible supports. 
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EQUATIONS OF MOTION 


Figure (^represents the single mass Jeffcott rotor mounted in damped 
elastic supports. In the Jeffcott model, the shaft is considered as a 
massless elastic member and the rotor mass is concentrated in a disc 
mounted at the center of the span. The shaft in term is supported in 
linear bearings which are mounted in damped flexible supports. 

Neglecting rotor acceleration and the disc gyroscopics, the governing 
equations of motion for the rotor, bearings, and support system in complex 
notation reduce to the following (See Appendix A for derivation): 

M 2 Z 2 + C Z 2 + C.Z e - QZ 2 + (K - IojC. )Z = M 2 e o) 2 e la>+ (I) 

z ^ s * is T z s is z p 

C.Z. - C.Z + K.Z, - (K - i ujC . ) Z = 0 (2) 

b j is b b s is 

M1Z1 + CiZi + K1Z1 - C.Z - (K - ! aiC . ) Z = 0 ( 3 ) 

1 ii ii j 5 s l S 


where 


Z = Z 2 - Z. - Zi 

s 2 j 1 


relative shaft deflection. 


If the internal damping C. and the aerodynamic cross coupling term Q 
are excluded from the above equations then the system will be stable (26). 

After the initial transient motion has damped out, it may be assumed 
that the system steady-state motion is circular synchronous precession. In 
this case the displacements are related to the velocity and acceleration 
vectors as follows: 


Z, = A . e ' U1+ 

Z. = iujZ. (4) 

Z. = luZ. = - w 2 Z. 

i i i 


* I I I ustrat i ons begin on page 38. 
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where A. is in general complex. 

The differential equations of motion may be reduced to a set of algebraic 
equations for the determination of the rotor steady-state motion. 


(K - 
s 

M 2 w 2 + iC u))Ao “ K A. - K Ai = M?e uj 2 
s * sj s 1 

(5) 

-K A 2 

+ (K, + K + IwCJA. + K Ai = 0 
b s b j s 1 

(6) 

-K A 2 
s *■ 

+ K A. + (Ki + K - Mia ) 2 + icoCi )Ai = 0 
s j 1 s 1 L 1 

(7) 


ROTOR AMPLIFICATION FACTOR 


Consider the steady-state orbit of the flexible rotor on rigid 
supports. The rotor amplitude Is a function of both the rotor and bearing 
stiffness and damping characteristics. Assuming Aj is zero, the relative 
journal bearing complex amplitude from Eq. 7 is given by 


K (K + 1C - i ojC. ) 

A. = -®_2 £ A 2 

J (K + K. ) 2 + (uC, ) 2 
s b b 

Solving Eq. 8 for the rotor amplitude yields 

( K 2 “ M 2 co 2 - i o>C ^ ) 

A 2 = M 2 e w 2 — 

U (K 2 - M 2 <o 2 ) 2 + ( <juC 2 ) 2 


( 8 ) 


(9) 


where 


K 2 

c 2 


K. K (K + K, ) + K (wC.) 
b s s b s b 

(K + K, ) 2 + (uC, ) 2 
$ b b 


K 2 C, 
s b 


(K, + K ) 2 + (ojC ,) 2 
b s b 


+ C 


2 


s 


The rotor displacement vector Z 2 may be expressed in terms of the 
absolute displacement R 2 and the phase angle <p as follows 




- R 2 e 


i (urf - <p) 


( 10 ) 
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where 


Ro - 


M 2 e u ur 


~\JTk~2 ~ M 2 <o 2 )^ + (coC 2 ) 
coC 2 


<t> = tan 1 


K 2 — M 2 oj 2 


The above results are similar to the rotor amplitude and phase angle 
results for the single mass flexible rotor on rigid supports as shown by 
Thomson (27). 

The rotor undamped, or natural critical speed is given by 



For the case of a lightly damped rotor system on rigid supports the 
maximum rotor amplitude will occur at approximately the rotor critical 
speed and the dimensionless rotor amplitude or amplification factor at 
the critical speed is given by 

R 2 k 2 

A = — | =— p- (12) 

e '<o=co co C 2 
u c c £ 


Examp I e I 

Consider a 97 lb. disc centered on a uniform massless elastic shaft 

K. 

as shown in Fig. (I). Assume that the bearing stiffness ■^ L is 500,000 
Ib/in and that the effective shaft stiffness at the disc station is 
333,000 Ib/in. Assuming light damping, the total stiffness K 2 is given by 

K 2 = — _ 1 x °- 555 x iq1 . 2 a 250,000 Ib/in 

K + K, (I + 0.333) 1 0 6 
s b 

The rotor critical speed is 

= 1,000 rad/sec 
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or N c = 9,550 RPM. 

If the rotor damping C g is assumed to be 15 lb - sec/ in and the bearing 
damping coefficient C^/2 is 80 ib — sec/in, then the effective system damp- 
ing coefficient C 2 is approximately given by 


C + 
s 


K 2 C. 
s b 


(K, + K ) 2 

D S 


|.5 + (0-333) 2 x IQ 12 x 160 
(I. 333) 2 x I0 12 


25 I b - sec/ in 


The ampl ification factor at the rotor critical speed is given by 
K 2 


CR 


A = 


250,000 _ 


u c C 2 1,000x25 


= 10.0 


The ampl ification factor of 10 represents a very I ightly damped rotor 
system and indicates that the rotor amplitude at the critical speed will be 
10 times the rotor unbalance eccentricity e . 


ROTOR RESPONSE ON DAMPED FLEXIBLE SUPPORTS 


Solution of Eq. 7 for the case of synchronous precession for the 


shaft relative deflection Z g yields 


Z s = CZ 2 - Z x ) 


K b (K b * V + { “ C b )2 + '“C b K s 


(K b + K s )2 + (<JjC b )2 


(13) 


Hence, in terms of the general coefficients C 2 and K, 


(Z ? - Zi ) 


K 2 + iu)(C 2 - C ) 


(14) 


The simultaneous equations for the absolute shaft and support housing 
motion reduce to the following 


□<! + K 2 - Mju) 2 + iw(C 2 + Cj - C s )]A 1 + C-iw(C 2 - C ) - K 2 ]A 2 = 0 

(15) 


6 



( 16 ) 


C-K 2 - i to (C2 “ C s )]Ai + [K2 _ M2U) 2 + I C0C2DA2 = 


If the damping terms are neglected, then the natural frequencies 
of the system may be determined by the expansion of the determinant of 
coefficients. The resulting frequency equation may be expressed as follows 


u 1>2 /“ c 



I + K 
2M 


± 



I + K 
2M 


+ 


1)2 _ i< 

2 M 


(17) 


where 



Figure 2 represents the dimensionless critical speeds vs. the dimen- 
sionless support stiffness factor K for various values of support to rotor 
mass ratios. Note that the incorporation of the flexible support with the 
rotor bearing system causes two critical speeds to occur; one which is 
higher and one which is lower than the original rotor critical on rigid 
supports. 

To solve for the complex support and rotor amplitudes A) and A 2 , Eq. 
16 and 17 may be expressed as follows: 


Ca u + lb ij > j = F r 


2 ; l 


(18) 


Multiplying Eq. 19 by the complex inverse matrix of coefficients and ex- 
panding yields 


F 1 ai2 + ibj2 

F 2 a 2 2 + ,b 22 


(19) 


Where 


A = d p + id. 

d r = ( K 2 - M 2 w 2 )(K 1 - Mju> 2 ) - K 2 M 2 W 2 — C 1 C 2 W 2 - o) 2 C s (C 2 - 
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dj = CjwO^ " M 2 w^ ) + C2 aJ (K i - Mjo) 2 - l^w^) + C co ( K 2 + N^u^) 


Expanding Eq. 20 

^l a 22 * ^2 a 12 + * t ^ 1 b 2 2 * ^2^12^ 


At = 


d + Id . 
r I 


( 20 ) 


In this case only an external unbalance excitation force F 2 is acting 
on the shaft and no external exciting force Fj is assumed to be present on 
the support system. For example, an excitation for Fj may be transmitted 
to the rotor system through the support structure by vibrations of auxiliary 
or adjacent equipment. 


F 2 Cai 2 d + b 12 d. + i(b 12 d - a 12 d.)J 
Ai i 1 ! — 


d 2 + d. 2 
r 1 


(21 ) 


Assume Aj Is of the form 


Aj A jr i Aj . 


( 22 ) 


The complex support amplitude Z\ after some complex algebraic manipulations 
Is given by 


Z 1 = ^e 


T cot 


Rje 


i ( cot - 


Si) 


( 23 ) 


where 


R 1 = V A i r 2 + A li 2 ' = +an “ 1 (j L ) 

r 


If the shaft damping coefficient C g is considered small in comparison 
to the effective damping coefficient C 2 than the system displacements and 
phase angles are given as follows 


R 1 = M 2 e u co 2 


I K 2 2 + (uC 2 ) 2 


dr 2 + d, 2 


( 24 ) 


and the phase angle of the support motion relative to the rotating unbal- 
ance load is given by 
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Bl = tan -1 




K 2 d p 4 coC 2 d . 


Since the complex rotor support motion Zi is given by 

i (cot - 0 1 ) 

Zi = Xj + I Y i = Rie 

Then for example, the horizontal and vertical components of 
motion are given by 


( Xj | / K 2 2 + (wC 2 ) 2 jcos(wt - Bj)) 

hi) =Me ^V d 2 7 d 2 (sinOot- M) 

v r i 

In a similar fashion, the complex rotor amplitude Z 2 is given by 


an + ibn . 


Z 2 - M 2 e u “ d 7 id, 


i cot 


r i 

After some manipulation, Eq. 28 reduces to the following 

/( K\ + K 2 • M^co 2 ) 2 + ( ( C x t C 2 )co) 2 1 (cot - 0 2 ) 


Z 2 = M 2 e u co 2 “ 


d 2 4 d . 2 
r i 


where 


0 2 = tan” 1 


( K] 4 K 2 — M^co 2 )d . - (Cj 4 C 2 )cod^ 
( K x 4 K 2 - M x o) 2 )d 4 (C x 4 C 2 )wd. 


The relative journal displacement is given by 


Z j ■ Z * - Z > - 


Where the relative shaft deflection is 
(Z 2 - Zj) 


K_ 


[K - K 2 - iu>C 2 ] 


Solving for the journal displacement 


(25) 


the support 


(26) 


(27) 


(28) 


(29) 


(30) 


(31 ) 
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i (cot - 8 . ) 

Z. = R. e J 

J b 

where, 

R. = Mo6 u 2 
J u 

II V, 1 S / 

(32) 

and the phase angle 3j between the journal amplitude and rotating un- 
balance force is given by 

(Kj - M 1 o) 2 )d. - coC i d \ ojCo i 

H 1 + +an-i lr — — (33) 

(Kj - M 1 o> 2 )d i + tuCid . j \s K2 j 

FORCES TRANSMITTED 




The magnitude of the resultant forces transmitted through the bearings 
and the support are of considerable interest to the designer from a stand- 
point of bearing life and system isolation. It is desirable to minimize 
the forces transmitted through the supporting structure and foundation so 
that other machines or piping systems are not excited. The magnitude of the 
force transmitted through the bearings is given by 



and the force transmitted through the support system is given by 


Fj = Rj y Kj 2 4 - (uC i ) 2 


An indication of the effectiveness of the support system in attenuating 
the forces transmitted to the foundation is the support dynamic transmissi- 
bllity factor TRD which will be defined as the ratio of the magnitude of the 
transmitted support force to the rotating unbalance load. If the dynamic 
transml ss I b 1 1 i ty is less than I, then the support system possesses good 
attenuation characteristics. Analysis has shown that if the support 
housing impedance character i st ics, which are determined by the housing mass. 
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stiffness and damping, are mismatched to the rotor-bearing system then, 
under certain speed conditions the dynamic transmi ssi bi I i ty may exceed I. 


The dynamic transmi ss i b i I i ty for the support is defined as 
F i /(K 2 2 + (uC 2 ) 2 )(K 1 2 + (C^) 2 ) 


TRD = 


(36) 


^ 2 e u ^ 


d 2 + d. 2 
r i 


If it is assumed that the rotor is operating well above any of the 
system critical speeds then the dynamic transmissibi I ity is approximately 
given by 


TRD = 



(K 2 2 + (wC 2 ) 2 )(K 1 2 + (C^) 2 ) 


M 2 2 M 2 2 


(37) 


The above expression leads to the well known conclusion that to 
minimize the forces transmitted through the support for supercr It ica I 
speed operation in the Jeffcott model, the support damping should be 
zero and the support stiffness should be as light as possible (28). 

This is a highly undesirable design practice for several reasons since 
large rotor amplitudes and forces transmitted may be encountered when pass- 
ing through the rotor critical speeds, and also the rotor system would be 
extremely shock sensitive and particularly susceptible to self-excited 
whirl Instability under such conditions. 

A compromise support damping coefficient should be selected to either 
minimize the rotor amplitudes or the forces transmitted over the operating 
speed range and also be sufficient to insure adequate rotor stability. 


ANALYSIS OF SYSTEM UNBALANCE RESPONSE - TUNED SYSTEM 


Figure 3 represents a computer generated plot of the dimensionless 
rotor relative amplitude versus the dimensionless rotor speed for the 
case of K = M = I. This relative rotor amplitude is equivalent to the 
motion monitored by a proximity probe mounted in the casing measuring 
the rotor motion at the center span. This system represents a tuned 
condition in which the support stiffness ratio K is equal to the support 



mass ratio M. With no support damping in the system, the tuned support 
will cause the relative rotor amplitude to be zero at a speed correspond- 
ing to the rotor critical speed with rigid supports. The Introduction 
of support mass and flexibility has caused two critical speeds to appear 
in the system; one above and one below the rigid support rotor critical. 

Note that when the support damping is relatively low the amplitudes at 
the two criticals becomes extremely high. 

As the dimensionless support damping ratio C increases from 0.01 to 
10 the rotor amplitudes at the system critical speeds decrease while the 
amplitude increases at a speed corresponding to the rigid support critical 
speed (u)/gj c = I). Note that in this case the damping value of 10 appears 
to be close to an optimum value for the minimization of the resonance 
amplitudes. If the support damping is further increased from 10 to 50, 

Fig. 5 indicates that there will be only one critical speed present in 
the system which will correspond to the rigid support critical. Although 
the damping of C = 50 is over 5 times the optimum value, the maximum 
amplitude is only 1/3 the rigid support value of 10. As the support damp- 
ing approaches infinity, the rotor amplitude will asymptotically approach 
10 . 

Figure 4 represents the absolute dimensionless rotor motion for various 
values of support damping ratio and is similar to Fig. 3. It should be 
noted that the damping coefficient of 10 also appears to be close to the 
optimum damping for the absolute motion as well as the relative motion. 

It is of interest to note that the various damping lines all intersect 
at a common point P in the plot of absolute as well as relative rotor motion. 
If the rotor amplification factor A is 100 (implying light rotor damping) 
then there will be two common points of intersection P and Q on the response 
plots (see Fig. 10) similar to that shown by Den Hartog for damped vibration 
absorber (15). The intersection points P and Q will occur at speeds res- 
pectively below and above the rigid support critical speed. The rotor 
amplitude may be minimized for the case of the absolute rotor motion by 
selecting the damping such that the slope of the response curve is zero 
at point P, and zero at point Q to minimize the rotor relative motion. 
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Figure 5 represents the phase angle between the rotating unbalance 
vector and the absolute rotor displacement vector for various damping 
coefficients. The phase angle for the single mass rotor on rigid supports 
(Jeffcott model) increases as the speed increases from 0 to 90 degrees 
at the critical speed and asymptotically approaches 180 degrees as the 
rotor speed greatly exceeds the critical speed. The phase angles of the 
rotor on damped flexible supports has a considerably different behavior 
from that of the rigid support rotor. For light values of supporting damping 
(C = 0.01), the phase angle increases rapidly to 180° as the system passes 
through the first critical speed and drops to almost 60° as it passes 
through the second critical speed. As the speed greatly exceeds the 
highest critical speed, the phase angle again approaches 180°. The phase 
angle of 180° indicates that the rotor mass center lies along the rotor 
spin axis. As the support damping coefficient is increased beyond 5 for 
the case of the tuned system, the reduction in phase angle above the first 
critical speed Is suppressed. This phenomena of phase angle reversal 
above the first critical speed has been observed experimentally. 

Figure 6 represents the support amplitude versus speed for various 
damping values and indicates that with very light support damping there 
will be large support resonances. As the damping is increased beyond 
C = 10 the resonances are suppressed and the amplitude is only slightly 
greater than I . For C = 50 there is only a smal I peak observed in the 
support system which occurs at a speed corresponding to the rigid support 
critical speed. The addition of high damping (C > 50) freezes the support 
and limits its motion drastically. 

Figure 7 represents the support phase angles versus speed ratio for 
various values of support damping. The phase angle for light damping 
(C = 0.01) is zero at low speeds and goes to 180 degrees as it passes 
through the first critical and then shifts to 330° upon passing through 
the second critical speed. If the rotor damping is light (A = 100) the 
support phase angle will approach 360° after passing through the second 
critical speed. Note that the various damping lines intercept at three 
points. The first node point represents the first system critical speed. 
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the second node point represents the rigid support critical speed and the 
third node point represents the second critical speed on flexible supports. 
In the discussion of the single mass flexible rotor presented in vibration 
texts (27) the phase change is only shown from zero to 180 degrees. In 
more complex systems with flexible supports, the phase change may vary 
between 0 and 360 degrees. For example In multimass systems the authors 
have observed phase changes of n times 180 degrees where n represents the 
number of system critical speeds. The measurement of rotor and support 
phase angles have been neglected and limited data has been reported in the 
literature. This is an extremely useful variable which when incorporated 
with displacement measurements can be used in balancing flexible rotors 
or impedance calculations of the support system. 

Figure 8 represents the dimensionless bearing forces transmitted for 
the tuned system. The dimensionless force transmitted is obtained by 
dividing by the transmitted force corresponding to the value at the crit- 
ical speed of the original rotor on rigid supports. Because of the light 
shaft damping the force transmitted curves are similar in appearance to 
the displacement curves. Note that for the support damping coefficient of 
C = 10 the forces transmitted to the bearings are only 17 percent of value 
transmitted for the rotor bearing system on rigid supports. 

Figure 9 represents the force transmitted through the bearing supports 
to the foundatation or base for various values of supporting damping. With 
a very lightly damped support system, (C = 0.01) the support amplitude and 
force transmitted will be particularly high at the first critical speed 
where the bearing and support motions are in phase. At the second critical 
speed, the support amplitude is lower than the amplitude attained at the 
first critical speed. This is because the bearings and support motions 
are out of phase which enables the bearing damping to help attenuate the 
support motion. It is of interest to note from Fig. 8, for the tuned 
rotor system, the bearing force transmitted at (to/ui^) = I with an undamped 
support system is zero. Figure 9 shows that the correspond i ng force trans- 
mitted through the support system at u/oj = I has been reduced to only I0£ 
of the rigid support value. 


14 



The force transmitted for an undamped support system at a speed 
ratio of four is approximately 10 % of the rigid support value. This con- 
dition would be desireable if it were possible to accelerate through the 
critical s, thereby avoiding the large steady-state amplitudes and forces 
developed. 

The near optimum damping of 10 increases the support forces transmitted 
in the supercritical speed region to ~50% of the rigid support value and 
the overdamped support system (C = 50) has increased to nearly 80%. Hence, 
the support damping introduced to suppress the system resonances will cause 
the forces transmitted to Increase in the supercritical speed region. 

If the system is designed to operate over the entire speed range shown, 
then the near optimum value of damping (i.e., C = 10) for suppressing the 
rotor absolute amplitude also produces the most desirable attentuation of 
forces to the system support structure. 

OPTIMUM DAMPING FOR TUNED SYSTEM 

From the observation of the computer generated displacement and force 
transmitted plots it is apparent that there exists an optimum damping to 
either minimize the rotor amplitudes or the forces transmitted over the 
entire speed range. 

For example to minimize the absolute rotor motion as shown in Fig. 4 
or the relative rotor motion shown in Fig. 3, the method of (16) may be 
used in which the damping is chosen so that the slope of the amplitude 
curve is zero at points P and Q respectively. In the tuned system where 

K/M = I for light rotor damping (A = 100), the rotor amplitudes at points 

P and 0 are independent of the support damping as shown in Fig. 10 and 

can be shown to be equal to 

X 2 = x 2 /e u l Pi p - ■'■HTIm < 58) 

Therefore with the tuned system illustrated with a mass ratio of M = I, 
the maximum amplitude at P or Q will be 1.732 times the rotor unbalance 
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eccentricity. The optimum damping may be selected so that the tangent to 
the amplitude curve at either point P or Q has a zero slope. By selecting 
the optimum damping in this fashion it is seen that the maximum amplitude 
in the system will not exceed the value given by Eq. 38. Thus it is readily 
apparent that to minimize the rotor response over a given speed range, the 
support mass should be kept as light as possible. 

After considerable algebraic manipulation (28) the optimum damping 
coefficient for both points P and Q is given by the following expression 

^2 _ 4M 3 ip 3 - 3M 2 (4M + 5 ?ip 2 + M( I 2 2 + I 3M + 8)0 - M(l + 2M) 2 ^ 

- 1 2Mip z + 8(1 + 2M)t jj - (I t 2M) 

where, 

c = Ci/C c = Cj/Cz x I/2A = Cx 2 ^- 

4i - 0 1 2 or ft 2 ^ depending on whether the value calculated is for point 

P or Q respecti vel y. 


and 


0,2 = /rT ^ T 

I + A t 2M 


o 2 2 = _2lZH_ 

. . . . /I + 2M - I 

For example, when M = I and for the first node, P: 

^ = Qj2 = — {S_ - 0.634 

i + AS 


and 


C 2 = 0.447 


Hence 
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Cl 

— = 0.688 for point P 

c opt 


In a similar fashion 


Cl, 

= 0.559 for point Q 

c 'opt 
Example 2 

As an example of the application of the tuned support design criteria 
consider the rotor of Example I mounted in flexibly supported bearing housings 
which weigh 48.5 lbs and have a stiffness of 125,000 Ib/in. The total 
support weight and stiffness Kj is given by 


Wi = 2 x 48.5 = 97 lb 

Kj = 2 x 125,000 = 250,000 Ib/in 


Hence, 


M = Mi/M 2 = I .0 
K = Kj/K 2 = 1.0 


The critical damping coefficient C 


C = 


2712 _ 500,000 Ib/in _ 

= - — tt~t - 500 


c oj c 1,000 rad/sec 


is given by 
I b-sec/i n 


Thus the support damping coefficients required to make the slope of 
the rotor amplitude curve zero at points P and 0 are respectively given 
as fol lows 


Ci I = 0.688 x C = 344 lb-sec/ in 
Ip c 

Ci I = 0.559 x C = 279.5 Ib-sec/in 
1 'q c 
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These calculations are valid only for the case of zero damping on 
the rotor and in the bearings (i.e., A = °») and only for the tuned system 
(i.e., K = M). For a more realistic solution, a value of A = 10 was chosen 
and numerous cases were then programmed on a digital computer to arrive at 
a value of optimum amplitude and required damping. This approach is dis- 
cussed in the next section of this paper but the results for the tuned 
system are very nearly the same as the results arrived at analytically for 
the case of A = » and are presented in Fig. II. 

The results shown in Fig. II are approximately correct for systems 
having moderate to light damping on the rotor (i.e., 10 < A < <*>). Note 
that the smaller the mass ratio M is, the lower the peak response will be 
and also the lower the required support damping will be. For example, if 
the mass ratio is 0.1, then the maximum dimensionless amplitude will be only 

I. l and the required damping ratio will be 5 as compared to a value of 13.6 
for an M ratio of I. Figure 12 is a response plot for the tuned system 

K = M = 0. I which i I I ustrates the val idity of the results plotted in Fig. 

II. The response curve for a damping ratio of 5 passes almost horizontal 
through the node point and has the low ampl itude ratio as indicated by 
Fig. II. 

Example 3 

Consider a rotor system simi lar to Example 2 in which the rotor rigid 
support amplification factor A = 10. 

For a tuned support system the dimension support damping coefficient 
is obtained from Fig. II for M = I as follows 

C = C 1 /C 2 = 13.6 

where C 2 is given as 25 lb-sec/ in (Example I). 

The re fore, 

C x = 13.6 x C 2 = 340 lb-sec/ in 
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Note that this value is approximately the same as the value given In Example 
2 for the required damping at point P corresonding to A^. 

This indicates that each support must have 170 lb-sec/ in damping to 
achieve the optimum response of about 1.7 times the unbalance level of the 
rotor. 

Next consider a tuned support with a mass and stiffness ratio of 0.10 
(see Fig. 12). Correspond i ng support weight and stiffness are given as 
fol lows 

Wj = 9.7 I b/ i n 

Kj = 25,000 Ib/in 

The required damping is thus found from Fig. II to be 
C = 5.0 


or 


Cj = 5 x 25 = 125 lb-sec/ in 

Thus only 62.5 lb-sec/ in damping per support is required to obtain an 
optimum response of l.l times the unbalance level of the rotor. This value 
of I . I is in comparison to a maximum response of 10 times the unbalance 
level for the rigidly mounted rotoi — bearing system. 

OPTIMIZATION OF SUPPORT DAMPING FOR OFF-TUNED CONDITIONS 

In general it is not possible or necessarily desirable to have a 
tuned support system. The support to rotor mass ratio is usually dic- 
tated by design cons i derat ions and can be varied only within certain 
ranges. Figure I I shows that for best reduction of rotor ampl itude, the 
support mass should remain as light as possible. However, it will be 
shown that even with high mass ratio support systems the rotor amplitudes 
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can be attenuated by a factor of 5 by proper selection of the stiffness 
and damping coefficients. 

To evaluate the optimum damping for off-tuned conditions the computer 
program was run for various support mass and stiffness ratios and each of 
these for various damping coefficients. For example. Fig. 13 represents 
the amplitudes at the rotor first and second critical speed for various 
mass ratios with a dimensionless stiffness ratio of K = I as the mass ratio 
and damping are varied. The solid lines represent the amplitude at the 
second critical speed and the dotted lines represent the amplitude at the 
first critical speed. With moderate support damping ratios it is observed 
that as the mass ratio increases the amplitude at the first critical reduces 
while the amplitude at the second critical increases. The optimum damping 
was selected as the intersection of the amplitudes at the first and second 
critical for a particular value of damping. For example, the lowest optimum 
amplitude point on the plot is given by a damping ratio of 10 and produces 
an amplitude ratio of about 1.5. Several plots similar to Fig. 13 were 
produced and the results were then crossplotted to obtain plots of amplitude 
versus damping ratio. 

Figure 14 represents the maximum rotor amplitude vs. support damping 
ratio for various values of dimensionless support stiffness for a rotor 
bearing system with a low support mass ratio of 0.01. Figure 14 shows 
that for this particular case, the lowest amplitude is achieved by a low 
support stiffness ratio of K = 0.01 which is of the same order as the mass 
ratio. With this low support stiffness, there is a wide range of support 
damping (re. C = I -*■ 6) that can be used to achieve the low level of 
rotor response. 

Thus, under proper design conditions the support damping may be 
allowed to vary by a considerable amount without Impairing the rotor 
performance. As the support stiffness ratio increases, the maximum rotor 
amplitude response also increases and the required support damping must 
be larger. For example, if the support stiffness ratio increases from 
0.01 to 2.0, the optimum damping required increases by a factor of 15 from 
approximately 2 to 30. 
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No+e also that for high stiffness ratio support systems, the per- 
missible range of the support damping coefficient is very narrow, and 
that either a reduction or an increase of damping beyond the optimum 
value will result in a rapid gain in rotor response. 

It is also of interest to note that if a high support stiffness (K = 2) 
is used in conjunction with a low value of support damping (C < 2) then the 
rotor response will be worse than the original rotor response on rigid 
supports (A = 10). 

Figure 15 represents the maximum rotor response vs. support damping 
for a high support mass ratio system (M = 2). It is obvious from the 
comparison of Figs. 14 and 15 that the high mass ratio support system is 
less desirable. The minimum rotor amplitude that can be achieved is X 2 /E = 2 
with a tuned support where K = M = 2 and a support damping coefficient of 
C = 20. (Also see Fig. II on the tuned system.) As the support stiffness 
ratio is reduced, the rotor response curve increases in the optimum damping 
reg i on . 

If it is not possible to incorporate a high value of support damping 
Into the system (C = 20), then the rotor amplitude can still be reduced 
to 40/f of the original rotor response by a low support damping value of 
C = I and a reduced support stiffness ratio of K = 0.7. For low values of 
support damping, if the support stiffness Increases beyond K = 0.7, the 
rotor response rapidly increases. 

A series of plots similar to Figs. 14 and 15 were produced for various 
mass ratios in order to determine the optimum rotor response for off-tuned 
support conditions. Figure 16 represents the rotor maximum amplitude vs. 
the support mass ratio for various values of support stiffness with optimum 
damping. 

For the case of A = 10, Fig. 16 illustrates that the lowest amplitude 
can be achived with a low mass ratio support system. With a high mass ratio 
support system such as M = 5, the rotor amplitude X 2 can be reduced from 10 
to 2.8 by means of a tuned support stiffness of K = 5.0 and optimum damping. 
Note that as the support stiffness becomes very light, the maximum rotor 
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amplitude increases to 7.5. 

At a low value of support mass (M = 0.1), the rotor amplitude increases 
as the support stiffness increases. The required optimum damping required 
with the low mass ratio support is found to be given approximately by the 
relationship that C = I 5K for M < 0.2. 

TRANSIENT ANALYSIS 

The previous discussion has been concerned only wi fh the steady-state 
response of the rotor due to unbalance and has not considered the rotor 
initial transient motion. As discussed previously, the damped flexible 
support system is impoitant, not only from the standpoint of reduction of 
synchronous unbalance response, but also in the control of self excited 
vibrations such as caused by internal friction, aerodynamic excitation, etc. 
Therefore to investigate the general rotor motion and also to provide a 
check on the steady-state analysis, the rotor equations of motion were 
integrated forward in time on the digital computer using a modified 4th 
order Runge-Kutta integration procedure. This procedure is of importance 
particularly if the analysis is extended from a linear bearing or support 
system to include a nonlinear hydrodynamic damper bearing as presented in 
Ref. 13. 

The dimensionless rotor and support transient orbits were automatically 
computer plotted by a Cal comp plotter with the following dimensionless 
parameters 

X = x/e , Y = y/e 
u' 1 u 

Figure 17 represents the initial transient orbit of a 96.6 lb rotor 
similar to Example I in highly damped support (C = 43) for the first 12 
cycles of shaft motion. The support mass ratio and the support stiffness 
ratio are both approximately the same (0.10) which represents a tuned 
system. Because of +he excessive support damping, the maximum force 
transmitted to the support is 2.16 times the unbalance force while the 
force transmitted to the bearings is reduced by about 40 %. The magnifications 
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of the force to the support would be highly undesirable for applications 
such as aircraft jet engines. For example, various investigators have 
observed that such a situation occurs with the hydrodydynamic squeeze film 
bearing when operating at excessive eccentricity ratios (29). 

Figure 18 represents the transient orbit for the same rotor system 
except that the support damping has been reduced by a factor of 100 from 
Ci = 1,000 Ib-sec/in.to 10 Ib-sec/in. In this case, the maximum force 
transmitted through the support is less than 16$ of the rotating unbalance 
force and the bearing force transmitted is below 9%. This orbit is 
analogous to a suddenly applied unbalance such as a blade loss in an engine. 
Although the forces transmitted have been greatly reduced with the low 
stiffness and damping support system, the rotor has developed a large 
initial transient motion of over 10 times the unbalance eccentricity and 
this transient motion is not readily damped out. 

In Fig. 19, the rotor transient motion is depicted with an optimum 
damping coefficient of C = 5.5 for minimum rotor response as determined 
fromlhe steady-state analysis. The transient response is rapidly suppressed 
after seven cycles of shaft motion to produce a small stable synchronous orbit. 
The transmitted forces to the bearings and support are nearly balanced 
to achieve approximately a 15% attenuation of the unbalance load. 

SUMMARY AND CONCLUSIONS 

The equations of motion for a single mass rotor-bearing system on 
damped flexible supports have been derived and studied considering both 
a steady-state and transient type analysis. Design charts for both tuned 
and off-tuned support conditions have been presented. 

The analysis may be summarized by the following general statements. 

I. The critical speed response of the single mass Jeffcott model 
rotor may be completely eliminated by means of a low mass ratio flexible 
support with optimum damping. In this case the rotor steady-state ampli- 
tude of motion over the entire speed range will only be slightly more than 
the rotor unbalance eccentricity. 
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2. The support mass ratio should be kept as light as possible 
to achieve minimum rotor amplitude. 

3. The rotor ampl itude may be considerably attenuated even for high 
mass ratio support systems by tuning the support stiffness such that 

K = M and i ncorporat i ng optimum damping for the tuned conditions. 

4. With a low mass ratio support system, the required value of 
optimum damping is not critical and can vary by a actor of 10 without 
appreciably effecting rotor performance. 

5. As the mass ratio increases the required value of optimum damping 
increases rapidly and the permissible range of vacation of the support 
damping diminishes. 

6. The off-tuned support (K M) can be designed to produce a con- 
siderable improvement in system response in comparison to the rotor on 
rigid supports. 

7. If insufficient damping is incorporated in the support then the 
resulting rotor steady-state amplitude may be larger than the original 
rotor response for support stiffness values K > I. 

8. If there is excessive support damping (C > 20) with a low mass 
ratio support (M = 0.1), then the forces transmitted through the support 
may exceed the unbalance forces (TRDS > 1.0). 

9. Although the steady-state analysis shows that the rotor amplitude 
will be small for an underdamped (C < 0.50) low mass ratio support system, 
the orbital analysis shows that a large initial transient motion can be 
generated due to the suddenly applied unbalance force and that this motion 
is not readily attenuated. 

10. The transient program indicates that a rotor with an underdamped 
and low stiffness ratio support system Is highly susceptable to self 
excited non-synchronous precession or whirling. 

11. The optimum damping based on minimization of the rotor steady- 
state ampl itude for both tuned and off-tuned conditions produces a satis- 
factory transient response from the standpoint of rapid reduction of the 
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initial transient motion, improved system stability and reduction of the 
forces transmitted. 
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APPENDIX A 


DISCUSSION 

A high speed rotor shaft may be considered as a continuous elastic 
member with variable mass and inertia properties along its length. The 
rotor shaft usual ly has attached to it such components as turbine or 
compressor blades, impeller disks, or spacer assemblies or seals. If 
the axial dimensions of each rotor component is small in comparison to the 
overal I length of the rotor, then each component may be treated as a con- 
centrated mass with a polar moment of inertia equivalent to that of the 
original component. If the mass of the components are large In comparison 
to the shaft mass connecting the components, then the shaft weight can 
be neglected or considered to be located at the mass stations. If the polar 
moment of Inertia of each section is ignored, then the stations may be 
considered as concentrated masses, rather than distributed In the plane 
of the rotor element. However, if the sections whirl In a plane, perpen- 
dicular to the spin axis then the gyroscopic moments do not act on the 
system and hence the equations reduce to the same as If point masses were 
assumed. 

The position vector of the nth mass center is given by 

\ + 7 J *1 +ir u 

where 

= vectoral bearing support deflection 

<5. = vectoral journal deflection 

J 

5 g = vectoral shaft deflection 

= displacement of mass center from the shaft centerl ine 
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The system being analyzed has been reduced to a single mass rotor 
mounted in idealized linear bearings and the bearings are in turn mounted 
on damped, elastic supports. By considering only small deflections, the 
spring rate of the flexible, massless rotor shaft may be considered to be 
I i near. 

The rotor disk (see Fig. I) is considered to whirl in a plane and 
hence no gyroscopic moments are acting on the system. The orthogonal support 
and bearing spring rates are assumed symmetric and no cross coupling terms 
are considered to be acting at the support housings. The aforementioned 
assumptions allow the equations of motion of the system to be written as 
total differential equations. 

DERIVATION OF EQUATIONS OF MOTION 

A. I Ki nematics 

The position vectors to the mass stations are given by 
m l : bearing M /0 

housing mass P = Xjn^ t Yjn^ (A. I ) 

m 2 : rotor mass 

M 2 /0 

P = (X 2 + e cos0)n + (Y 2 + e sin9)n (A. 2) 

z u x *■ u y 

The velocities of the mass stations are given by 
_ Mj/O . 

V = XirT + Yirf (A. 3) 

1 x 1 y 

M 2 /O • • * • 

- (Xo - e 0sin0)"n + (Yo + e 0cos0)rT (A. 4 ) 

z u x z u y 
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A. 2 Kinetic Energy 


The kinetic energy of the system is given by 


4 i 


i = 1 


(M. V . 

i i 


V 


1 3 3 

+ T I 1 • -W-U) . 

2 i=l j=l 'J 1 J 


(A. 5) 


neglecting the gyroscopic couples acting on the disc, the system kinetic 
energy reduces to 


T - Ml ( X i 2 t Yi 2 ) t ~ M 2 [(X 2 - e u 0sin0) 2 + (Y 2 - e u 0cos0) 2 U 


+ i V 2 


(A. 6) 


A. 3 Potential Energy 

The potential energy of the system is composed of the sums of the 
potential energy of the flexible shaft, the potential energy of the bearings, 
and the energy of support structure as fol lows 

V = j CK s (X 5 2 t Y s 2 ) + K b (X. 2 t Yj 2 ) t KjfX! 2 t Y! 2 )J (A. 7) 


where 



A. 4 Dissipative Energy 

The system dissipative energy consists of the damping functions provided 
by the bearing support system, the bearings, and the external and internal 
rotor damping and the aerodynamic rotor cross coupling as described by 
Alford (4). 

D = i (Cj (Xj 2 t Y! 2 ) t C.(X. 2 t Y. 2 ) + C (X 2 2 t Y 2 2 ) 

2 b j j s z z 

• • * » • * 

+ C.[X 2 t Y 2 t 2 uj (Y X - X Y )] t Q(Y 2 X 2 - X 2 Y 2 )} (A. 8) 

IS S 5 S S S 
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The internal damping function is dependent upon the rotor precession rate 
and can cause self excited whirl instability when the rotor is operated 
above the critical speed (26). Alford has demonstrated that the aerodynamic 
cross coupling stiffness term can also cause rotor Instability when the 
rotor speed is supercritical . When the system dissipation function is com- 
prised of only the first three terms, the system Is inherently stable. 

A. 5 Lagranges Equations 

The governing equations of motion are obtained from Lagranges Equations, 
which state: 


c— : - 

3q 


— + - D 
3q r 9q r 



where 


(A. 9) 


L = T - V 

The total number of equations of motion obtained will be equal to the 
number of degrees of freedom of the system which is seven and are given as 
fol lows. 


Rotor 


X 2 : M 2 X 2 + C s X 2 + C. (X 2 - X! - X fc ) + « s (X 2 - X 1 - X ) 

+ QY 2 + o)C.(Y 2 - Yi - Y^) = M^Ja^costiot) + asin(oot)) 


(A. 10) 


Y 2 : M 2 Y 2 + C Y 2 + C.(Y 2 - Yi - Y.) + K (Y 2 - Y x - Y ) - QX 2 

si J ^ J 


(A. I I ) 


- toC. (X 2 - Xi - Xj ) = M^ u (co 2 s i n (tot) - asin(wt)) 
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The equations A. 10 to A. 15 may be vectorially combined by representing 
the displacements in complex notation as follows 



Bear i n 



Ul 


(A. 20) 
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NOMENCLATURE 


A 

cr 

A j 

Ai 

A 2 

C 

C b 

c 

c 

C. 

i 

c 

s 

c l 

C 2 

e 

u 

Fi 

F b 

K 

K b 

K 

s 

Kl 

K 2 

M 

Mj 

m 2 

P 

0 

Q 


Amplification factor at rigid support critical 

Complex bearing amplitude, in 

Complex support amplitude, in 

Complex rotor amplitude. In 

Damping ratio = Cj/C 2 (DIM) 

Bearing damping, lb-sec/ in 

Critical damp i ng coeff ic lent, Ib-sec/in 

Rotor internal damping, Ib-sec/in 

Absolute shaft damping, Ib-sec/in 

Support damping, Ib-sec/in 

Effective rotor-bear i ng damping, Ib-sec/in 

Rotor mass eccentricity, in 

Force transmitted to foundation, lb 

Force transmitted to bearing housing, lb 

Stiffness ratio, Ki/K 2 

Bearing stiffness, lb/ in 

Rotor-shaft stiffness, lb/ in 

Support stiffness, I b/ i n 

Effective rotor-bear ing stiffness, lb/ in 

Mass ratio, = Mj/M 2 (DIM) 

Support mass, lb-sec 2 / in 
Rotor mass, ib-sec 2 /in 


K 2 

—7 r- (DIM) 
oj C 2 


1st node point on response plot 


Rotor cross-coupling stiffness, I b/ i n 


2nd node on response plots 
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r 2 

T 

TRD 

V 

V 

Wi 

X 

X 

s 

Xi 

x 2 

X. 

J 

V 

s 

Yi 

Y 2 


Zi 

z 2 

a 

Pi 

$2 

6 b 

e 


Y 

€ 

<{> 


Rotor absolu displacement amplitude 
Kinetic energy 

Transmi ssi bi I ity = F j/ (M 2 e u oo 2 ) 

Potential Energy 
Velocity, in/sec 
Support weight, I b 
Def i ned as X 2 

Shaft relative displacement in x - direction 

Support displacement in x - direction 

Rotor absolute displacement in x - direction 

Journal relative displacement in x - direction 

Shaft relative displacement in y - direction 

Support displacement in y - direction 

Rotor absolute displacement in y - direction 

Journal relative displacement in y - direction 

Complex shaft relative amplitude 

Complex support amp I itude 

Complex rotor amplitude 

Rotor angular acceleration, rad/sec 2 

Phase angle of support motion relative to rotor unbalance, DEG 
Phase angle of rotor motion relative to rotor unbalance, DEG 
Phase angle of bearing motion relative to rotor unbalance, DEG 
Angular displacement, rad 
Defined as K/M (DIM) 

Damping ratio = cj/c (DIM) 

Rotor absolute amplitude phase angle, DEG 
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Moment of inertia 


$ 

X Optimum amplitude for tuned system 

¥ Defined as fij 2 or n 2 2 when calculating required damping at point 

P or p respect i vel y 

u Rotor angular velocity, rad/sec 

ui >2 Rotor system critical speeds, rad/sec 

u c Rigid support critical speed, rad/sec 

^ 1 * &2 Speeds at which the node point P and Q occur on response plots 
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DAMPER 


Figure I 



Schematic Diagram of Single Mass Rotor on Damped Elastic Supports 
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ROTOR CRITICAL SPEEDS w. 9 /<u r (DIM.) 



RELATIVE AMP. RATIO CX2-X1 )/E [DIM. 3 

0.000 1.000 2.000 3.000 4.000 5.000 6.000 



FREQUENCY RATIO (W/WC) 

Figure 3. Dimensionless Relative Rotor Amplitude Vs Speed Ratio for Various Values 
of Support Damping for a Tuned Support System, K = M - | 
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Figure 4. Absolute Rotor Motion with a Tuned Support System for Various Values 
of Support Damping, K = M = I, A = 10 
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FREQUENCY RATIO CW/WC) 









AMPLITUOE RATIO < 



Figure 6. Support Amplitude Vs Speed for Various Values of Support Damping 
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Figure 7. Phase Angle of Support Motion Relative to Rotor Unbalance for Various 
Values of Support Damping 





BERRINGS VS FREQUENCY RATIO 



45 


FREQUENCY RATIO CW/WC) 





FORCE TRRNS. / FTRoWC FOR RIGID SUPPORT 
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Figure 9. Dimensionless Force Transmitted to Foundation Vs Speed Ratio for 
Various Values of Support Damping 






AMPLITUDE RATIO [X 2 /e] DIMENSIONLESS 



Figure H. Optimum Support Damping and Maximum Rotor Amplitude Vs Mass Ratio for 
A ~ 10 
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Figure 12. Rotor Amplitude Vs Speed for a Low Mass Ratio Tuned Support System 
for Various Values of Support Damping K = M = 0.1, A = 10 
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FOR VARIOUS STIFFNESS RATIOS 
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Figure 14. Rotor Maximum Amplitude Vs Damping Ratio for Various Values of Stiffness 
Ratios for a Low Mass Ratio Support, M = 0.01, A = 10 
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FOR VARIOUS STIFFNESS RATIOS, K 



Figure 15. Rotor Maximum Amplitude Vs Oamping Ratio for Various Values of Stiffness 
Ratios for a High Mass Ratio Support, M = 2, A = 10 
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WITH OPTIMUM SUPPORT DAMPING 
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Figure 16. Rotor Maximum Amplitude for Various Values of Stiffness and Mass 
Ratio with Optimum Support Damping 
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ABSOLUTE ROTOR MOTOION 


N = 30000 RPM M = 0. 100 

K = 0.092 C = 43.916 

W2 = 96.60 LB. KB « 500,000 LB/IN 

KS » 500,000 LB/ IN CB = 100.0 LB-SEC/ IN 

DC = 0.5 LB-SEC/ IN W1 * 9.66 LB. 


CD = 0.5 LB-SEC/IN K1 = 25,000 LB/IN 
QRC = 0.5 LB/IN Cl * 1000.0 LB-SEC/IN 
TRDB = 0.585 AND OCCURS AT 0.58 CYCLES 



X-DIR. 

Figure 17. Dimensionless Transient Motion of an Unbalanced Rotor for 12 Cycles on 
Overdamped Supports (K = M = 0. I , C = 44 ) 
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Figure 18. 
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X-DIR. 

Dimensionless Transient Rotor Motion with Underdamped Flexible Supports 
for 12 Cycles (K = M = 0,10, C = 0.44) 
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ABSOLUTE ROTOR MOTOION 
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Figure 19. Dimensionless Rotor Motion with Optimum Steady-State Damping Showing 

the Steady-State Orbit After Seven Cycles of Running Speed (K = M = 0.1, 
C = 5.5) 
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